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BREAKING SPACES AND FORMS FOR THE DPG METHOD AND 
APPLICATIONS INCLUDING MAXWELL EQUATIONS 


C. CARSTENSEN, L. DEMKOWICZ, AND J. GOPALAKRISHNAN 

Abstract. Discontinuous Petrov Galerkin (DPG) methods are made easily implementable 
using “broken” test spaces, i.e., spaces of functions with no continuity constraints across 
mesh element interfaces. Broken spaces derivable from a standard exact sequence of first 
order (unbroken) Sobolev spaces are of particular interest. A characterization of interface 
spaces that connect the broken spaces to their unbroken counterparts is provided. Stability 
of certain formulations using the broken spaces can be derived from the stability of analogues 
that use unbroken spaces. This technique is used to provide a complete error analysis of DPG 
methods for Maxwell equations with perfect electric boundary conditions. The technique 
also permits considerable simplifications of previous analyses of DPG methods for other 
equations. Reliability and efficiency estimates for an error indicator also follow. Einally, 
the equivalence of stability for various formulations of the same Maxwell problem is proved, 
including the strong form, the ultraweak form, and various forms in between. 


1. Introduction 

When a domain Q is partitioned into elements, a function in a Sobolev space like i7(curl, Q) 
or i7(div, Q) has continuity constraints across element interfaces, e.g, the former has tangen¬ 
tial continuity, while the latter has continuity of its normal component. If these continuity 
constraints are removed from the space, then we obtain “broken” Sobolev spaces. Discontinu¬ 
ous Petrov Galerkin (DPG) methods introduced in [15, 17] used spaces of such discontinuous 
functions in broken Sobolev spaces to localize certain computations. The studies in this paper 
begin by clarifying this process of breaking Sobolev spaces. This process, sometimes called 
hybridization, has been well studied within a discrete setting. For instance, the hybridized 
Raviart-Thomas method [5, 32] is obtained by discretizing a variational formulation and then 
removing the continuity constraints of the discrete space, i.e., by discretizing first and then 
hybridizing. In contrast, in this paper, we identify methods obtained by hybridizing first and 
then discretizing, a setting more natural for DPG methods. We then take this idea further 
by connecting the stability of formulations with broken spaces and unbroken spaces, leading 
to the first convergence proof of a DPG method for Maxwell equations. 

The next section (Section 2) is devoted to a study of the interface spaces that arise when 
breaking Sobolev spaces. These infinite-dimensional interface spaces can be used to connect 
the broken and the unbroken spaces. The main result of Section 2, contained in Theorem 2.3, 
makes this connection precise and provides an elementary characterization (by duality) of 
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the natural norms on these interface spaces. This theorem can be viewed as a generalization 
of a similar result in [33]. 

Having discussed breaking spaces, we proceed to break variational formulations in Sec¬ 
tion 3. The motivation for the theory in that section is that some variational formulations 
set in broken spaces have another closely related variational formulation set in their unbroken 
counterpart. This is the case with all the formulations on which the DPG method is based. 
The main observation of Section 3 is a simple result (Theorem 3.1) which in its abstract form 
seems to be already known in other studies [25]. In the DPG context, it provides sufficient 
conditions under which stability of broken forms follow from stability of their unbroken rela¬ 
tives. As a consequence of this observation, we are able to simplify many previous analyses 
of DPG methods. The content of Sections 2 and 3 can be understood without reference to 
the DPG method. 

A quick introduction to the DPG method is given in Section 4, where known conditions 
needed for a priori and a posteriori error analysis are also presented. One of the conditions is 
the existence of a Fortin operator. Anticipating the needs of the Maxwell application, we then 
present, in Section 5, a sequence of Fortin operators for i7(curl, iF) and i7(div,iF), 

all on a single tetrahedral mesh element K. They are constructed to satisfy certain moment 
conditions required for analysis of DPG methods. They fit into a commuting diagram that 
helps us prove the required norm estimates (see Theorem 5.1). 

The time-harmonic Maxwell equations within a cavity are considered afterward in Sec¬ 
tion 6. Focusing first on a simple DPG method for Maxwell equation, called the primal DPG 
method, we provide a complete analysis using the tools developed in the previous section. To 
understand one of the novelties here, recall that the wellposedness of the Maxwell equations 
is guaranteed as soon as the excitation frequency of the harmonic wave is different from a 
cavity resonance. However, this wellposedness is not directly inherited by most standard 
discretizations, which are often known to be stable solely in an asymptotic regime [29]. The 
discrete spaces used must be sufficiently fine before one can even guarantee solvability of the 
discrete system, not to mention error guarantees. Furthermore, the analysis of the standard 
finite element method does not clarify how fine the mesh needs to be to ensure that the stable 
regime is reached. In contrast, the DPG schemes, having inherited their stability from the 
exact equations, are stable no matter how coarse the mesh is. This advantage is striking 
when attempting robust adaptive meshing strategies. 

Another focus of Section 6 is the understanding of a proliferation of formulations for 
the Maxwell boundary value problem. One may decide to treat individual equations of 
the Maxwell system differently, e.g., one equation may be imposed strongly, while another 
may be imposed weakly via integration by parts. Mixed methods make a particular choice, 
while primal methods make a different choice. We will show (see Theorem 6.3) that the 
stability of one formulation implies the stability of five others. The proof is an interesting 
application of the closed range theorem. However, when the DPG methodology is applied 
to discretize these formulations, the numerical results reported in Section 7, show that the 
various methods do exhibit differences. This is because the functional settings are different for 
different formulations, i.e., convergence to the solution occurs in different norms. Section 7 
also provides results from numerical investigations on issues where the theory is currently 
silent. 
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2. Breaking Sobolev spaces 

In this section, we discuss precisely what we mean by breaking Sobolev spaces using a mesh. 
We will define broken spaces and interface spaces and prove a duality result that clarifies the 
interplay between these spaces. We work with infinite-dimensional (but mesh-dependent) 
spaces on an open bounded domain f] with Lipschitz boundary. The mesh, denoted by 
is a disjoint partitioning of into open elements K such that the union of their closures 
is the closure of f7. The collection of element boundaries dK for all K e is denoted by 
dQh- We assume that each element boundary dK is Lipschitz. The shape of the elements is 
otherwise arbitrary for now. 

We focus on the most commonly occurring first order Sobolev spaces of real or complex¬ 
valued functions, namely iL(div, f7), and iL(curl, f7). Their broken versions are de¬ 

fined, respectively, by 

= {ne L\n) : u\k e H\K), K e Qu) =0 H\K), 

i/(curl, Qh) = {Ee : E\k e i/(curl,iL), K e Qh] =0 

KeQh 

H{diY,Qh) = W e : G\KeH{diY,K),KEQh] = H H{diY,K). 

KeQh 

As these broken spaces contain functions with no continuity requirements at element inter¬ 
faces, their discretization is easier than that of globally conforming spaces. 

To recover the original Sobolev spaces from these broken spaces, we need traces and inter¬ 
face variables. First, let us consider these traces on each element K in 

ueH'^{K), 

tr^rl.T E = i'^K X E) X nK\dK E e if (curl, K), 

tr^j .1 E = nK X E\sk E e if (curl, if), 

= a\sK ■ riK (7 6 if(div,if). 

Here and throughout uk denotes the unit outward normal on dK and is often simply written 
as n. Both uk and these traces are well defined almost everywhere on diF, thanks to our 
assumption that dK is Lipschitz. The operators trgrad? f^curi,T? fi^curiH? perform the 

above trace operation element by element on each of the broken spaces we defined previously, 
thus giving rise to linear maps 

trg,^d : H\Qh) - ft tre„ri,T : H{cni\,Qh) - ft ^“'/'(curl, dff), 

KeQh KeQh 

trcurl,H : ^f(curl,iffe) ^ P[ if"^/^(div, dff), traiv : ff(div, ^ P[ H~^^^{dK). 

KeQh KeQh 

It is well known that these maps are continuous and surjective (for the standard definitions 
of the above codomain Sobolev spaces, see e.g., [29]). An element of YlKeiih H~^^‘^{dK) is 
expressed using notations like n • d or (even when d itself has not been assigned any 
separate meaning) that are evocative of their dependence on the interface normals. Similarly, 



4 


C. CARSTENSEN, L. DEMKOWICZ, AND J. GOPALAKRISHNAN 


the elements of the other trace map codomains are expressed using notations like 

nxE = E^e P[ dK), {nxE)xn = Eye P[ dK). 

KeQh KeQh 

Next, we need spaces of interface functions. We use the above trace operators to define 
them, after cautiously noting two issues that can arise on an interface piece / = dK^ n dK~ 
shared by two mesh elements K- in First, functions in the range of trgrad when restricted 
to F is generally multivalued, and we would like our interface functions to be single valued in 
some sense. Second, the range of the remaining trace operators consists of functionals whose 
restrictions to / are in general undefined. The following definitions circumvent these issues. 

= trgrad Qh) = ticud.T -^^(curl, J7), 

dfih) = t^curiH -ff(curl, fl), = ti'div -ff(cliv, fl). 

If Qh consists of a single element, then H^/‘^{dQh) equals but in general tig^ad £ 

tVgj-ad (and similar remarks apply for the other spaces). We norm each of the above 

interface spaces by these quotient norms: 


(2.1a) 

(2.1b) 

(2.1c) 

(2.1d) 


ll-^T||_H--i/2(div,ai7h) 

-®Hllif-i/2(curl,ai2ft) 

Il^n||i/-i/2(5i7^) 


inf Ikllnqr?), 

inf , l|-E^llff(curl,J7)) 

£;eH(curl,«)ntr-ll ^{Bh} 

inf „ ll-®lli?(curl,r?)) 

EeHicvir\,Q)r.tr-^^^-^{ET} 

inf _ Ikllff(div,i7)- 

creif (div,i7) ntr^.^ {dn } 


These are indeed quotient norms because the infimums are over cosets generated by kernels 
of the trace maps. E.g., if w is any function in such that trgrad'^ = then the 

set where the minimization is carried out in (2.1a), namely n equals the 

coset w + YlKeOh where H^{K) = ker(tr^g^^). Note that every element of this coset 

is an extension of u. For this reason, such norms are also known as the “minimum energy 
extension” norms. For an alternate way to characterize the interface spaces, see [34]. 


Remark 2.1. The quotient norm in (2.Id) appeared in the literature as early as [33]. The 
word “hybrid” that appears in their title was used to refer to situations where, to quote [33], 
“the constraint of interelement continuity has been removed at the expense of introducing 
a Lagrange multiplier.” The quote also summarizes the discussion of this section well. The 
above definitions of our four interface spaces are thus generalizations of a definition in [33] 
and each can be interpreted as an appropriate space of Lagrange multipliers. 


We now show by elementary arguments that the quotient norms on the two pairs of trace 
spaces 

{H^/^{dK),H-^/^{dK)} and {H-^/^{cutI, dK), dK)}, 

are dual to each other. The duality pairing in any Hilbert space X, namely the action of a 
linear or conjugate linear (antilinear) functional x' e X' on x e X is denoted by (x'^x)x'xx 
and we omit the subscript in this notation when no confusion can arise. We also adopt the 
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convention that when taking supremum over vector spaces (such as in the next result) the 
zero element is omitted tacitly. 

Lemma 2.2. The following identities hold for any n-a in and any u in 


(2.2a) 

(2.2b) 




llif-V 2 (^K) - sup 




Mn^i^idK) = sup 


|<n ■ a, u)\ 


= sup 
ueH^/^(dK) 

sup 


\(^n,u)\ 


<7eH{div,K) lkllif(div,K) aneH-'^^/'^idK) W^nWH-y^idK) 

The next two identities hold for any E_^ in H~^/^{diY^dK) and any Ft in H~^/^{cml,dK). 

\(E^,F}\ 


(2.2c) 


\\H{cur\,K) 

\<E^Tt}\ 


(2.2d) 


\\F-t\\H-'^/2(div,eK) — sup 

FeH{cml,K) 

= sup 
FjeH-y^ (curl, 5i<r) II -^TII H-F2 (curl,5i^) 

IP II _ \(nxE,FT)\ 

I-frII h- 1/2(curl,5i<r) “ SUp 

EeH (curl,K) l|-C'llif(curl,/r) 

I<^h,A>I 


= sup - 

E-ieH-F^{div,dK) 


WE. 


Hllif-i/2(div,aK) 


Proof. The first identity is proved using an equivalence between a Dirichlet and a Neu¬ 
mann problem. The Dirichlet problem is the problem of finding a e H{div,K), given 
dn e H~^/‘^{dK), such that 


(2.3) 


(2.4) 


n ' a = n ' 

on dK^ 

— grad(div cr) + cr = 0, 

in K. 

w E H^{K) satisfying 


dw 


on 

on dK^ 

—div(gradrc) w = 0^ 

in K. 


It is immediate that problems (2.4) and (2.3) are equivalent in the sense that w solves (2.4) 
if and only if a = gradu; solves (2.3) and moreover ||rc||ifi(/c) = ||<7||i/(div,i^)- It is also obvious 
from the calculus of variations that among all H{div, iir)-extensions of (T„, the solution of (2.3) 
has the minimal H{div,K) norm (i.e., a is the “minimum energy extension” referred to 
earlier), so 


h|lif-i/2(ai^) — lk||ij(div,/C) - ll^llfi'i(A') 

l<^n,^^)| 


= sup 


:eH^(K) iTWH'^iK) 


sup 

vem{K) 


Kgradui, gradu)^ + {w,v)k 
\\vWm{K) 
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where we used the variational form of (2.4) in the last step. (Here and throughout, we use 
(•, •)k to denote the inner product in Lp‘{K) or its Cartesian products.) This proves the first 
equality of (2.2a). 

Next, analogous to (2.3) and (2.4), we set up another pair of Dirichlet and Neumann 
problems. The first problem is to find u in given any u g such that 


(2.5) 


u = on dK^ 

—div(gradR) u = 0^ in K. 


The second is to find r in i7(div,i^) such that 


( 2 . 6 ) 


( div r = on dK^ 

\ — grad(div r) + r = 0, in K. 


The solution u of (2.5) has the minimal H^{K) norm among all extensions of u into H^{K)^ 
i.e., \MHV2^eK) = MhHk)‘ Thus *>/||R||ifi/ 2 (aK) = <d^,R>/||i/||^i(K), so taking the 
supremum over all u in we obtain 


sup 

uEH^/‘^{dK) 


|<dn,d)| 

IIII if 1/2 (aX) 


^ sup 
uEm{K) 


|<dn,^)| 

ll'^llifi(A:) 


Since the reverse inequality is obvious from the definition of the quotient norm in the denom¬ 
inator, we have established the second identity of (2.2a). To prove (2.2b), we begin, as above, 
by observing that r is the solution to the Neumann problem (2.6) if and only if r = divr 
solves the Dirichlet problem (2.5). Moreover, ||T||if(div,A:) ^ \\'^\\h^{k)- Hence 


II'^11 if 1/2(dAT) “ ll'^lliii(A:) “ 


\^\\H{div,K) - sup 

pEH{div,K) 

= sup , 

pEH{div,K) \\p\\H(div,K) 


|(divT,divp)K + {r,p)K\ 

\\p\\H{div,K) 


where we have used the variational form of (2.6) in the last step. The proof of (2.2b) can 
now be completed as before. 

We follow exactly the same reasoning for the i7(curl) case, summarized as follows: On one 
hand, the norm of an interface function equals the norm of a minimum energy extension, while 
on the other hand, it equals the norm of the inverse of a Riesz map applied to a functional 
generated by the interface function. The minimum energy extension that yields the interface 
norm ||^H||if-i/2(div dK) solution of the Dirichlet problem of finding E e iJ(curl, K) 

satisfying 

{ nx E = on dK, 

\ curl(curl£') E = 0, in K, 


while the inverse of the Riesz map applied to the functional generated by E^ is obtained by 
solving the Neumann problem 


(2.8) 


n X (curlF) = E^, on dK 

curl(curlF) E = 0, in K. 
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Again, the two problems are equivalent in the sense that F solves (2.8) if and only ii E = 
curlF solves (2.7). Moreover, \\E\\h(cut\,k) = ll^llff(curi,i^)- Hence 


WE. 


H llH-l/2(div,5i<r) “ W^Wh{cut\,K) - II-^11 if (curl,iC) 


sup 

GeH{cut\,K) 


I (curl F, curl G)k + (^, G) 


K 


sup 

GeH{cut\,K) 


IIGII if (curl,K) 


IIII if (curl,K) 


The proof of (2.2c) follows from this. The proof of (2.2d) is similar and is left to the 
reader. □ 


Let us return to the product spaces like iL(curl, and iL(div, f2h). Any Hilbert 

space V that is the Cartesian product of various Hilbert spaces V{K) is normed in the 
standard fashion, 

v= Yi v{K), iii;|i7 = y; wvkWI^k^, 

KeQh KeQh 

where vk denotes the F-component of any v in V. The dual space V' is the Cartesian product 
of component duals V{Ky. Writing an ^ g W as l{v) = YjKeQh where Ik ^ ^(^)^ 

it is elementary to prove that ||^||y/ = YjKeQh ll^^lly(K)'’ 


(2.9) 


(sup 
\veV 


\^{v)\ 

Wv\\v 



sup 

VKeV(K) 


\G{vk)\ 

ll^is:lk(K) 


2 


Some of our interface spaces have such functionals, e.g., the function in H gives 

rise to £{v) = where 


(Gn^v)h — 2 '^^)if-i/2(^X)xif 1/2(aX) 5 


KE^]h 


is a functional acting on g which is the sum of component functionals = 

(du, ^’)if-i/2(^i^)xiiV2(aK) acting on vk = v\k over every K g Other functionals like 

are defined similarly. We are now ready to state a few basic relationships between 
the interface and broken spaces. As usual, we define H^{Q) = {i’ g H^{Q) : v\dQ = 0}, 
iL(div, C) = {r G iL(div, Q) : = 0}, and iL(curl, Q) = {F e iL(curl, Q) : = 0}. 
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Theorem 2.3. The following identities hold for any interface space function dn in H , 

u in E_^ in H~^l‘^{dYV,dK), and Fj- in H~^^‘^{cw:l,dK). 


(2.10a) 

(2.10b) 

(2.10c) 

(2.10d) 


Wr 


IIh-i/ 2(5|2^) — sup 


\(^n,u}h\ 




uem{QH) 

\(n ■ a, u)h\ 


P II II ’ 

creHidiv,^^) IPm(div,J7;,) 

IF II KEh.Gi.I 

FEH{cuT\,Qh) 11^ II if (curl,i?/,) 

\(ri X 


IIF- 


T||ii-i/2(curl,^i7/,) 


£;eif (cud,17/,) Il^l|if(curl,l7/,) 


For any broken space function v g t g iJ(div, f2h)^ cind F g iJ(curl, iljf), 


(2.11a) 

veH\Q) 

(o-n, V)h = 0 


(2.11b) 

T e ii(div, 17) 

! 

II 

o 

'^iie F/^{dQh), 

(2.11c) 

F 6 .^(curl, 17) 

^ <Eh,FX = 0 

6 ii-i/2(div, 17/,) 


Proof. The identities immediately follow from Lemma 2.2 and (2.9). The proofs of the three 
equivalences in (2.11) are similar, so we will only detail the last one. If F is in iL(curl, 12), 
then choosing any E g iL(curl, 12) such that trcuri,T E = E^ and integrating by parts over 
entire 12, 


(curlL;,F)r2 + {E,c\itIE)q = 0, 

because of the boundary conditions on E on 512. Now, if the left hand side is integrated by 
parts again, this time element by element, then we find that {E^^Ey^ = 0. 

Conversely, given that (E^^Ey^ = 0 for any E in iL(curl, 12/^), consider curlF g (T>(12)^)'. 
As a distribution, curlF acts on 0 g T>(12)^, and satisfies 

(curlF)((/)) = (F, curl 0)^2 = (curlF, 0)/, - (n x (jy^Ey^ = (curlF, 0)/,, 

where we have integrated by parts element by element and denoted 

(•, ■)h = L (■’ 

KeQh 

This notation also serves to emphasize that the term curl E appearing on the right-hand side 
above is a derivative taken piecewise, element by element. Clearly curl F|x is in L?{KY for all 
K E f2h since E g iL(curl, 12/^), so the distribution curlF is in L^(12). Having established that 
F 6 if (curl, n), we may now integrate by parts to get <n x F, -F>j/-i/2(div,5J7)xi/-i/2(curi,ar2) = 
(curlF,F)r2 + (F,curlF)r2 = (n x E^Ey^ = 0 for all E e iL(curl, 12). This shows that the 
trace {n x E) x n\er 2 = 0, i.e, E g iL(curl, 12). □ 

Remark 2.4. While and UKeO.H^^HSK) are dual to each other, our 

interface spaces ^md are not dual to each other in general. 
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Remark 2.5. Equivalences analogous to (2.11) hold with interface subspaces 

= trgrad Qh) = ticurin ^(curl, Q), 

H~^^^{cux\,dQh) = t^curi.T-ff(curl, i7), H~^^‘^{df2h) = tr^iv-ff(div, i7). 

By a minor modification of the arguments in the proof in Theorem 2.3, we can prove that 
for any v e g iJ(div, 17/^), and F g i7(curl, 


(2.12a) 

veH^{f2) <= 

(^n,v}h = 0 

< 

q> 

(D 

1 

(2.12b) 

r G iL(div, f7) (r • n, u}h = 0 

\/ue H^/^{dQh), 

(2.12c) 

FeH(cuil,n) <= 

X 

II 

o 



3. Breaking variational forms 

The goal in this section is to investigate in what sense a variational formulation can be 
reformulated using broken spaces without losing stability. We will describe the main result 
in an abstract setting first and close the section with simple examples that use the results of 
the previous section. 

Let Xq and Y denote two Hilbert spaces and let Yq be a closed subspace of Y. For 
definiteness, we assume that all our spaces in this section are over C (but our results hold also 
for spaces over M). In the examples we have in mind, Y will be a broken space, while Yq will be 
its unbroken analogue (but no such assumption is needed to understand the upcoming results 
abstractly). The abstract setting involves a continuous sesquilinear form bo : Xq x Y ^ C 
satisfying the following assumption. 

Assumption 1. There is a positive constant cq such that 

II II ^ I^o(t,2/)| w ^ ^ 

cofIIao ^ sup ^- Mx G Aq. 

yevo \\y\\y 

It is a well-known result of Babuska and Necas [1, 30] that Assumption 1 together with 
triviality of 

(3.1) Zo = [y ^Yo ho{x^ y) = ^ foi* sill x g Xq} 

guarantees wellposedness of the following variational problem: Given I g Yq' (the space of 
conjugate linear functionals on Yq), find x g Xq satisfying 

(3.2) ho{x, y) = £{y) Vy e Yq. 

When Zq is non-trivial, we can still obtain existence of a solution x provided the load func¬ 
tional i satisfies the compatibility condition i{z) = 0 for all z g Zq. In (3.2), the trial space 
Xq need not be the same as the test space Yq. 

To describe a “broken” version of (3.2), we need another Hilbert space X, together with a 
continuous sesquilinear form b : X x Y C. In applications Y and X will usually be set to 
a broken Sobolev space and an interface space, respectively. Define 

K (x, x), y) = bo{x, y) + b{x, y). 

Clearly b : X x Y ^ C is continuous, where 

(3.3) X = Xo X 1 



10 


C. CARSTENSEN, L. DEMKOWICZ, AND J. GOPALAKRISHNAN 


is a Hilbert space under the Cartesian product norm. Now consider the following new broken 
variational formulation: Given i find x g Xq and x e X satisfying 

(3.4) b{{x,x), y) = e{y) ^yeY. 

The close relationship between problems (3.4) and (3.2) is readily revealed under the following 
assumption. 

Assumption 2. The spaces To? ^5 and X satisfy 

(3.5) To = G T : h{x, y) = {) for all x e X} 
and there is a positive constant c such that 

(3.6) c||x|U ^ VxeX 

yeY My 

Under this assumption, we present a simple result which shows that the broken form (3.4) 
inherits stability from the original unbroken form (3.2). A very similar such abstract result 
was formulated and proved in [25, Appendix A] and used for other applications. Our proof 
is simple, unsurprising, and uses the same type of arguments from the early days of mixed 
methods [5, p. 40]: stability of a larger system can be obtained in a triangular fashion by first 
restricting to a smaller subspace and obtaining stability there, followed by a backsubstitution- 
like step. Below, ||5o|| denotes the smallest number C for which the inequality \bo{x^y)\ ^ 
C'llxjlxo IIz/IIe holds for all x e Xq and all y eY. 


Theorem 3.1. Assumptions 1 and 2 imply 


where ci is defined by 


Ci\\{x,x)\\x ^ sup 

yeY 


1 _ 1 1 

A ^0 ^ 


\b{ (x,x), y)\ 


IT 


Co 


Moreover, if Z = {y eY : b{{x,x),y) = 0 for all x g Xq and x in X}, then 

Z = Zq. 

Consequently, if Zq = {0}, then (3.4) is uniquely solvable and moreover the solution eompo- 
nent x from (3.4) eoineides with the solution of (3.2). 


Proof. We need to bound \\x\\xo and \\x\\x- First, 

II II ^ \bo{x,y)\ 

coWAxo ^ sup ——Jj- 

yoeYo My 

\bo{x,y) + b{x,y)\ 


by Assumption 1, 


^ sup 

yoeYo 


IT 


by Assumption 2, (3.5) 
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Next, to bound ||x||^, using (3.6) of Assumption 2, 


a PI s sup - sup 

yeY ||y||y yeY ||y||y 


^ ||^>o|| Ikllxo + sup 
yeY 


\b{ {x,x), y)\ 


l|y 


Using the already proved bound for \\x\\xo in the last inequality and combining, 


{x,x)fx = Mxo + 11 ®"^ 


yco^eV My J \c\Co J yeY My J 

from which the inequality of the theorem follows. 

Finally, to prove that Z = Zq, using (3.5), 

y E Z 5o(x, y) = 0 for all x e Xq and 5(x, y) = 0 for all x e X 
5o(x, y) = 0 for all x g Xq and y g Yq? 

which holds if and only if ^ G Zq. 


□ 


Remark 3.2. Note that in the proof of the inf-sup condition, we did not fully use (3.5). We 
only needed Yq ^ {y e Y : h{x^y) = 0 for all x g X}. The reverse inclusion was needed to 
conclude that Z = Z^. 


Remark 3.3. It is natural to ask, in the same spirit as Theorem 3.1, if the numerical solutions 
of DPG methods using discretizations of the broken formulations coincide with those of 
discretizations of the original unbroken formulation. A result addressing this question is 
given in [3, Theorem 2.6]. 

In the remainder of this section, we illustrate how to apply this theorem on some examples. 

Example 3.4 (Primal DPG formulation). Suppose / g and u satisfies 

(3.7a) —Xu = f in f7, 

(3.7b) u = 0 on df2. 

The standard variational formulation for this problem, finds u in such that 

(3.8) (gradix,grad^)r 2 = 

This form is obtained by multiplying (3.7a) by 'L’ g and integrating by parts over the 

entire domain f7. If on the other hand, we multiply (3.7a) by a 'z; g and integrate by 

parts element by element, then we obtain another variational formulation proposed in [18]: 
Solve for u in as well as a separate unknown g (representing the 

fluxes —n • grad u along mesh interfaces) satisfying 

(grad n, grad + <(T„,y>/, = {f,y)Q yyeH^{fih)- 


(3.9) 
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Figure 1 . Chains of implications of inf-sup conditions 

We can view this as the broken version of (3.8) by setting 

Xo = H\Q), Yo = H\Q), 

X = H-^l\dQh), Y = H\Qh), 

^o(M,y) = (grad M, grad y)/,, 6((T„,y) = (an,y)h- 

For these settings, the conditions required to apply Theorem 3.1 are verified as follows. 
Coercivity of 6 o (*5 *) on Yq Assumption 1 holds. 

Theorem 2.3, (2.10a) (3-6) of Assumption 2 holds with c = 1. 

Theorem 2.3, (2.11a) (3-5) of Assumption 2 holds. 

Noting that Zq = {0}, an application of Theorem 3.1 implies that problem (3.9) is wellposed. 
This wellposedness result also shows that (3.9) is uniquely solvable with a more general 
right-hand side / in 

An alternate (and longer) proof of this wellposedness result can be found in [18]. The 
classical work of [33] also uses the spaces and but proceeds to develop a 

Bubnov-Galerkin hybrid formulation different from the Petrov-Galerkin formulation (3.9). /// 

Example 3.5 (Many formulations of an elliptic problem). Considering a model problem in¬ 
volving diffusion, convection, and reaction terms, we now show how to analyze, all at once, 
its various variational formulations. The diffusion coefficient a = a~^ : C ^ is a sym¬ 

metric matrix function which is uniformly bounded and positive definite on C, the convection 
coefficient is /3 g L°^(f7)^ which satisfies div(a/3) = 0 , and reaction is incorporated through a 
non-negative 7 g L^{Q). The classical form of the equations on Q are a = agradR-ha/3R-h/i 
and — diver-hyR = /2 (for some given fi g L^(C)^ and /2 g L?{Q)) together with the bound¬ 
ary condition u\^q = 0. This can be written in operator form using 

aa — grad u 
div a — (3 • a — 

We begin with the formulation closest to the classical form. 

Strong form: Let x = (cr, u) be a group variable. Set spaces by 

Xo = iL(div, C) X iL^C), Y = Yo = x l2(C), 

and consider the problem of finding x g Xq, given f ^ Y^ satisfying Ax = /. We can 
trivially fit this into our variational framework (3.2) by setting bo to 

bi{x,y) = {Ax,y)a. 


(3.10) 


A 


a 


u 



aa — grad u — 
div a — 'yu 
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Unlike the remaining formulations below, there is no need to discuss a broken version of 
the above strong form as the test space already admits discontinuous functions. The next 
formulation is often derived directly from a second order equation obtained by eliminating a 
from the strong form. 

Primal form: First, set spaces by 

Xo = H\f2), Yo = H\f2), 

X = Y = 

Then, with b set to b^{an, v) = (n ■ a, vyh and 6o set to 

bQ{u,v) = {agT:adu,gradv)h+iaPu,gi&dv)h + {'^u,v)q, 

the standard primal formulation is (3.2) and its broken version is (3.4). 

Next, consider the formulation derived by multiplying each equation in the strong form by a 
test function and integrating both equations by parts, i.e., both equations are imposed weakly. 
It was previously studied in [9, 16], but we can now simplify its analysis considerably using 
Theorem 3.1. 

Ultraweak form: Set group variables x = y = (r,i;), x = (d^,ii), and 

Xo = X Yo = H{diY, Q) x 

X = HY\dQh) X H-Y\dnh), Y = H{dW, Hh) x H\f2h), 
bo{x,y) = {x,A*y)Q^, b^{x,y) = (an,v)h-(u,n ■ T)h. 

Formulations (3.2) and (3.4) with bo = bg and b = W are of the ultraweak type. 

The fourth formulation, well-known as the mixed form [5], is derived by weakly imposing (via 
integration by parts) the first equation of the strong form, but strongly imposing the second 
equation. 

Dual Mixed form: Set the spaces by 

Xo = i/(div, Q) X l2(17), To = i/(div, Q) x 

X = y = i?(div, Qh) X L^{Q). 

The well-known mixed formulation is then (3.2) with ^ 

{t,v)) = {aax)h + ('^,divr)/, - {^ux)h 

-h (diver,'?;)/, - {ju,v)h. 

Its broken version is (3.4) with b set to b^{u^ (u'^)) ^ —{u^n • r)/,. 

Note that the well-known discrete hybrid mixed method [5, 11] is also derived from b^. That 
method however works with a Bubnov-Galerkin formulation obtained by breaking both the 
trial and the test i7(div, f2) components, while above we have broken only the test space. The 
last formulation in this example reverses the roles by weakly imposing the second equation 
of the strong form and strongly imposing the first equation: 
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Mixed form: Set 

Xo = X Yo = x 

X = Y = L^{Qf X 

The dual mixed formulation is (3.2) with 6o = ^ 

h^{ (cr,ix), (t,p) ) = (aax)a - (gTadux)n - {l3ux)o 
- (cr,gradp)i7 - (7R,p)r2. 

and its broken version is (3.4) with b set to 6^(d^, (r, p)) = {an^v}h- 

The variational problem (3.2) with is sometimes called [4] the primal mixed form to 
differentiate it with the dual mixed form given by b^. The broken formulation (3.4) with 
b^ and b^ was called the mild weak DPG formulation in [6]. Their analysis can also be 
simplified now using Theorem 3.1. 

In order to apply Theorem 3.1 to all these formulations, we need to verify Assumption 1. 
This can be done for all the formulations at once, because the six implications displayed 
in Figure 1 are proved in [13] for the model problem of this example (thus making the five 
statements in Figure 1 equivalent). We will not detail this proof here because we provide 
full proofs of similar implications for Maxwell equations in Section 6 (and this example is 
simpler than the Maxwell case). To apply these implications for the current example, we 
pick a formulation for which Assumption 1 is easy to prove: That the primal form is coercive 
b^ u) ^ C'IIR follows immediately by integration by parts and the Poincare inequality 

(under the simplifying assumptions we placed on the coefficients). This verifies Assumption 1 
for the primal form, which in turn verifies it for all the formulations by the above chain 
of equivalences. Assumption 2 can be immediately verified for all the formulations using 
either (2.11) or (2.12). Together with the easily verified triviality of Zq in each case, we have 
proven the wellposedness of all the formulations above, including the broken ones. /// 

4. The DPG method 

In this section, we quickly introduce the DPG method, indicate why the broken spaces are 
needed for practical reasons within the DPG method, and recall known abstract conditions 
under which an error analysis can be conducted. 

Let X and Y be Hilbert spaces and let 6 : X x T ^ C be a continuous sesquilinear form. 
In the applications we have in mind, X will always be of the form (3.3) (but we need not 
assume it for the theory in this section). The variational problem is to find x in X, given 
i eY'^ satisfying 

(4.1) b{x,y)=i{y) \/yEY. 

The DPG method uses finite-dimensional subspaces ez X and Y^ ei Y. The test space 
used in the method is a subspace ^ Y^ of approximately optimal test funetions computed 
for any arbitrarily given trial space X^- It is defined by = Th{Xh) where : X^ Y^ 
is given by 

(4.2) 


{ThZ,y)Y = b{z,y) My e Yh. 



15 


Here (•,•)>" is the inner product in y, hence by Riesz representation theorem on the 
operator is well defined. The discrete problem posed by the DPG discretization is to find 
Xh e Xh satisfying 

(4.3) bixh,y)=e{y) Vyeryl 

For practical implementation purposes, it is important to note that can be easily and 
inexpensively computed via (4.2) provided the space is a subspace of a broken space. 
Then (4.2) becomes a series of small decoupled problems on each element. 

For a posteriori error estimation, we use an estimator fj that actually works for any in 
X/i, computed as follows. (Note that Xh need not equal the solution Xh of (4.3).) First we 
solve for Sh in by 

(4.4) {£h,y)Y = ^{y) -b{xh,y), ^y^Yh. 

Again, this amounts to a local computation if Y^ is a subspace of a broken space. Then, set 

V = Il4||y. 

When Yh is a broken space, the element-wise norms of Sh serve as good error estimators [20]. 
The notations r] and Sh (without tilde) refer to similarly computed quantities with x^ in place 
of Xh> An analysis of errors and error estimators of the DPG method can be conducted using 
the following assumption introduced in [26]. In accordance with the traditions in the theory 
of mixed methods [5], we will call the operator 77 in the assumption a Fortin operator. 

Assumption 3. There is a continuous linear operator 11 : Y Y^ such that for all e X^ 
and all v eY^ 

b{wh^ V — nv) = 0 . 

Theorem 4.1. Suppose Assumption 3 holds. Assume also that there is a positive eonstant 
Cl sueh that 

(4.5) Cl IIXIIX ^ sup \/xeX^ 

yeY \\y\\Y 

and the set Z = {y e Y : b{x^y) = 0 for all x e X} equals {0}. Then the DPG method (4.3) 
is uniquely solvable for x^ and the a priori error estimate 

(4.6) ||(^ — 3;/i||x ^ ^ inf — Zh\x (quasi-optimality) 

Cl 

holds, where x is the unique exaet solution of (4.1). Moreover, we have the following inequal¬ 
ities for any x^ in X^ and its eorresponding error estimator fj, with the data-approximation 
error osc(7) = ||7 o (1 — 7T)||y/. 

(4.7a) ci||x — Xh\x ^ ll^ll V + osc(7), (reliability) 

(4.7b) fj ^ ||6|| \\x - XhWx, (efficiency) 

(4.7c) osc(7) ^ ||6|| ||1 — iT|| min ||a; — z/^llx. 
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Here ||7T|| and ||6|| are any constants that satisfy ||^^||y ^ ||7T||||^||y and \b{w^y)\ ^ 
||6||||R;||x||i/||y, respectively, for all w e X and y e Y. To apply the theorem to specific 
examples of DPG methods, we must verify (4.5). This will usually be done by appealing to 
Theorem 3.1 and verifying Assumptions 1 and 2. The previous sections provided tools for 
verifying Assumptions 1 and 2. In the next section, we will provide some tools to verify the 
remaining major condition in the theorem, namely Assumption 3. 

Remark 4.2. A proof of Theorem 4.1 is available in existing literature. The a priori error 
bound (4.6) was proved in [26]. The inequalities of (4.7), useful for a posteriori error esti¬ 
mation, were proved in [7]. In particular, a reliability estimate slightly different from (4.7a) 
(with worse constants) was proved in [7], but the same ideas yield (4.7a) easily (for example, 
cf. [8, proof of Lemma 3.6]). 

Remark 4.3. The operator is an approximation to an idealized trial-to-test operator T : 
X ^ Y given by 

{Tx,y)Y = h{x,y) My e Y. 

li B \ X ^ Y' Is the operator defined by the form satisfying {Bx){y) = 6(x, y) for all x e X 
and y E Y^ then clearly T = Ry^B^ where Ry : T ^ T' is the Riesz map defined by 
{RYy){v) = {y^v)Y‘ In some examples [16], it is possible to analytically compute T and then 
one may substitute with the exaetly optimal test spaee = T{Xh). 

Remark 4.4. The above-mentioned trial-to-test operator T = Ry^B should not be confused 
with another trial-to-test operator S = {B')~^Rx of [2] (also cf. [22]): 

X Y' 


Rx 


Ry 


X' 


B' 


Y 


Application of S requires the inversion of the dual operator B'. 


5. Fortin operators 

The Fortin operator 77 appearing in Assumption 3 is problem specific since it depends 
on the form b and the spaces. However, there are a few Fortin operators that have proved 
widely useful for analyzing DPG methods, including one for Y = another for 

Y = 77(div, f7/i), both given in [26]. In this section, we complete this collection by adding 
another operator for Y = 77(curl, intimately connected to the other two operators. Its 
utility will be clear in a subsequent section. 

Since the Fortin operators for DPG methods are to be defined on broken Sobolev spaces, 
their construction can be done focusing solely on one element. We will now assume that 
the mesh is a geometrically conforming finite element mesh of tetrahedral elements. Let 
Pp{D) denote the set of polynomials of degree at most p on a domain D and let Np{D) = 
Pp-i{D)^ X X Pp_i(D)^ denote the Nedelec [31] space. For domains D a M^, n = 2,3, let 
Rp{D) = Pp-i{D)'^ -h xPp-i{D) denote the Raviart-Thomas [32] space. We use Bp to denote 
the L? orthogonal projection onto Pp{K). From now on, let us use C to denote a generic 
constant independent of hx = diamX. Its value at different occurrences may differ and may 
possibly depend on the shape regularity of K and the polynomial degree p. 
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Theorem 5.1. On any tetrahedron K, there are operators 

: H\K) - Pp+3{K), 
-^p+3 : H{cml,K) ^ Np+3iK), 
n^%:H{dw,K)^Rp+3{K), 

such that the norm estimates 


(5.1a) 

(5.1b) 

(5.1c) 

hold, the diagram 
0-^ H^{K)/R 


n\ 


grad 

P+3 


\^pT3^\\h(cui\,K) ^ C '|| F ||^( curl , K )) 

II -^p+3 ^l|j/(div,iC) ^ C'||g||j^(div,A') 


Pp+3iK)/R 


grad 


grad 


i?(curl, K) 


curl 


H{div, K) 


div 


nitz 


/ydiv 

^^p+3 


Np+3{K) 


L\K) 

np+2 


^ Rp+3iK) Pp+2{K) 


commutes, and these identities hold for any v e H^{K), E e i?(curl, K), and r e H{div, K): 


(5.2a) 

(5.2b) 

(5.2c) 

(5.2d) 

(5.2e) 

(5.2f) 


(q, n^+3 v-v)k = 0 
(n ■ a, = 0 

{q,n;r^E-E)K = o 

<(n X F) X n, n X (iT^"lj F - F)> = 0 
(q, iTp^s r - r)^ = 0 
<n • {n^+3 r - r), /i> = 0 


V(?ePp_i(F), 

M n-a e Pp+i(F:), 
Vg6Pp(F)3, 

V (n X F) X n e ip+i(i^) 

V(?eFp+i(F)3, 

V M e tr^ad ^P+2(-fs^)- 


Note that the duality pairings above must be taken in the appropriate spaces, as in (2.2). 
To provide a constructive proof of Theorem 5.1, we will exhibit Fortin operators. We will 
use the exact sequence properties of the finite element spaces appearing as codomains of the 
operators in the theorem. We cannot use the canonical interpolation operators in these finite 
element spaces because they do not satisfy (5.1). Hence we will restrict the codomains of 
our operators to the following subspaces whose construction is motivated by zeroing out the 
unbounded degrees of freedom. 


= {i; G Pp+ 3 (iF) : v vanishes on all edges and vertices of iF}, 
= 1+6 Np+ 3 {K) : t • F = 0 on all edges of K, 




f (pn- curl F = 0 Vd 6 F°;4(dF) and 

JdK 

I ((n X E) X n) ‘ r = 0 \fr e Ri(dK) i, 

JdK j 

|r 6 Rp+ 3 {K) : j^^4>n-T = 0 V(/> e P^+^idK)^ 
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Here t denotes a tangent vector along the underlying edge, Ri{dK) = {r : r\f e Ri{f) for all 
faces / of dK}^ and and defined as follows. To simplify notation, 

let Pp{dK) = tr^^ (the space of functions on dK that are polynomials of degree 

at most p on each face of dK) and let Pp{dK) = tT^^^^Pp{K). Let Pp'^{dK) denote the 
L^(di^)-orthogonal complement of Pp{dK) + P{){dK) in Pp{dK) and let Pp{dK) denote the 
L^(di^)-orthogonal complement of Pp{dK) in Pp{dK). The following result is proved in [26, 
Lemma 3.2]. 


Lemma 5.2. For any v e H^{K), there is a unique in satisfying 

(5.3a) (g, TTq V — v)k = 0 for all q e Pp-i(K)^ 

(5.3b) (n • (j, TTq p — p) = 0 for all n • a e tr^^ Rp^i{K)^ and 

(5.3c) II TT| v\\l‘^(k) + ^kW gradTT| p||l 2 (r:) ^ C (||p||l 2 (r:) + hxW gradp||^ 2 (^)) . 

We define TT^+^^ as a minor modification of the analogous operator in [26]. Given any 
V E TT^(Tf), first compute its mean value on the boundary 

then split v = vq where vq = v — has zero mean trace, and finally define 

(5.4) 77^:;.^ p = nl VO + meK{v). 

Lemma 5.3. TT^^^^p satisfies (5.2a)-(5.2b) and (5.1a) for any v e H^{K) 

Proof. Since TT^^^^p — v = TTq pq “ equations (5.3a) and (5.3b) immediately yield (5.2a) 
and (5.2b). To prove the norm estimate (5.1a), note that standard scaling arguments imply 

|T^| 

(5.5) II^^a:('^)IIl2(r:) ^ ^ ^i\\^\\‘L‘^{K) + ^1:11 Si*adp||^2(^)), 

(5.6) ||p - meK{v)\\LHK) ^ ChK\\ gradp||^2(^). 
for all V in H^{K). Combining (5.5) and (5.3c), we get 

II ^pTs ^\\k{k) ^ C (||p||^2(^) + hxW gradp||^2(^)) , 

while combining (5.6) and (5.3c), 

hxW grad(i7®+f ^^)IIl2(k) = ^kW gradilf (^; - meK{v))\\L^K), by (5.4) 

^ C (||^; - msK(.v)\\i^2(^K) + hK\\ grad(^; - msK(.v))\\i^2(^K)) by (5.3c) 

^ ChK\\ gr&dvWL^K) by (5.6). 

These estimates together prove (5.1a). □ 

The next lemma is proved in [26, Lemma 3.3]. It defines TT^^g exactly as in [26]. 


Lemma 5.4. Any a e Bp^^^{K) satisfying 
(5.7a) {q.(T)K = 0 

(5.7b) (n • (7, /i) = 0 


yqePj,+i{Kf, 

VAietr^ad-Pp+2(^), 
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vanishes. Moreover, for any r e H{diY,K), there is a unique funetion in 

satisfying (5.2e)-(5.2f). It also satisfies (5.1c). 

The remaining operator ^^+3 will be defined after the next result. It is modeled after the 
previous two lemmas, but requires considerably more work. 

Lemma 5.5. Any E e satisfying 

(5.8a) (^^E)k = 0 y^ePp{Kf 

(5.8b) (/a,n X E) = 0 ^iie tr^^tr Pp+i{Kf^ 

vanishes. 


Proof. Integrating by parts twice and using (5.8b), we have 


(5.9) 


fjn ‘ curl E = (curl E, grad fj)K = (ji x E, grad fj) = {) \/fj e Pp +2 {K). 

JdK 


In addition, by Stokes theorem applied to one face / of K, we have 


(5.10) 


I A^n-curl£’= | nE‘t = 0 \fnePo(dK), 
Jf Jdf 


since t = 0 on all edges by the definition of B^^l{K). The definition of also gives 


(5.11) 


f 

JdK 


) n • curl E = 0 V<^ e {dK). 


Since Pp+ 2 {SK) = +Pp_,_2(5i^)+ Po(5-f^)) equations (5.9),(5.10) and (5.11) together 

imply 


(5.12) 


n • curl E = 0 

^dK 


e Pp+ 2 {dK). 


Since n • curlE’ e Pp^ 2 {,SK), we thus find that n • curl£^ = 0 on dK. This implies that the 
tangential component of E on dK, namely E^ = {n x E) x n, has vanishing surface curl, so 
it must equal a surface gradient, i.e., E^ = grad^'z; for some ^ ^ Pp+Z (dK). Moreover, since 
Ej vanishes on all edges, v may be chosen to be of the form v = bfVp for some Vp e Pp{dK), 
where 5/ is the product of all barycentric coordinates of K that do not vanish a.e. on /. 

To use the remaining (as yet unused) condition in the definition of B'^l^K), note that the 
tangential component of the coordinate vector x, namely is in Ri{dK). Combining this 
with (5.8b), we find that for all ji e tr^^^ Pp+i(i^)^ and any g M, 


(5.13) 


0 = (n X £', /i) + Ej • X = £'t • (/i X n + nx^ 

JdK JdK 

= 5] f grad^(5;^^).(/ixn + A.x^) 

/ 

= divT(/i X n + ^cXt), 


where the sums run over all faces / of dK. For any fi e tr^^^ ^ Pp+i(iF)^, the function ji x n 
is in the Raviart-Thomas space on the closed manifold dK denoted by Kfj^i{dK). (Note that 



20 


C. CARSTENSEN, L. DEMKOWICZ, AND J. GOPALAKRISHNAN 


unlike Ri{dK)^ this space consists of functions with the appropriate compatibility conditions 
across edges of dK.) The surface divergence map 

divT : Rl+i{dK) ^\wE Pp{dK) : f w = 0 

I JdK 

is surjective. Hence the term divT(/i x n + kxt) appearing in (5.13) spans all of Pp{dK) as 
ji and are varied. Choosing fi and k, so that divT(/i x n + kxj) = Vp^ we conclude that Vp 
vanishes and hence Ej = grad^(5jPp) = 0, i.e., 

(5.14) n X E = 0 on dK. 


Next, setting (j) = curlr in (5.8a) and integrating by parts, we obtain 

(5.15) r r*curl£^ = 0 Mr e Ppj^i{K)^. 

JK 

From (5.12) and (5.15), it follows that r = curlF^ is in and furthermore, r satis¬ 

fies (5.7a) and (5.7b). Hence, by Lemma 5.4, r vanishes. Thus curlF' = 0 and consequently 
E = gradp for some v e Pp^^{K). Furthermore, by (5.14), we may choose v = bxVp-i for 
some Vp-i G Pp-i{K)^ where bx is the product of all barycentric coordinates of K. Then 
(5.8a) implies 

gYad{bKVp-i) ■ (f) = \ bKVp-i ■ div(f) = 0 \f(p e Pp{K)^. 

JK JK 

It now follows from the surjectivity of div : Pp{K)^ Pp-i{K) that p^-i, and in turn 
E = grad(5KPp_i), vanishes on K. □ 

The next lemma defines the operator It will be useful to observe now that for any 

(5.16) I ((n X E) X n) ‘ r = 0 Vr e Ri(dK) | ((n x E) x n) • Xj = 0. 

JdK JdK 

Indeed, while the forward implication is obvious, the converse follows from (5.10). This 
shows that the condition that appears both in the definition of and in (5.16) above, 

actually amounts to just one constraint. 

Lemma 5.6. Given any E e i7(curl, iF), there is a unique B^^^E in Bp^^{K) satisfy¬ 
ing (5.2c)-(5.2d). 

Proof. We need to estimate n = dim Bp^^{K). First, note that since Po{dK) n is 

a one-dimensional space of constant functions on diF, 

dimP°:^^{dK) = dimPp+2{dK) - dim(P;+2(^^) + Po{dK)) 

= dimPp+2{dK) - dimP^+2iSK) - dimPo{dK) + 1 

= 6p 11. 


The tangential component E -t of any E e Np^s{K) is a polynomial of degree at most p -h 2 
on each edge, so E • t represents p -h 3 constraints per edge. Hence, counting the number of 
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constraints in the definition of 

n = dimBp"3(itr) ^ dim A^p+3(itr) — 6{p + 3) — dim(P°jj_'^(5itr)) — 1 
= dim A/p+3(K) — Q{p + 3) — {Qp + 11) — 1 
where we have used (5.16). Thus, 

(5.17) n ^ dim A^p+3(7f) — 12p — 30. 

Next, we count the number of equations in (5.2c)-(5.2d), namely 

m = dim(tr^^i^^ Pp^i{Kf) + dim Pp{Kf 
= 2dimPp^i{dK) — 6{p + 2) + dim Pp{K)^ 

= dimNp^s{K) - Q{p + 2) - 6{p + 3). 

This together with (5.17) implies that m = dimNp^^i^) — 12p — 30 ^ n Thus, the sys¬ 
tem (5.2c)-(5.2d), after using a basis, is an m x n matrix system of the form Ax = d, 
where x g is the vector of coefficients in a basis expansion of 77^+3 E and d is the 
right-hand side vector made using the given E. By Lemma 5.5, nul(A) = {0}. Hence 
m ^ n = rank(A) -h nul(74) = rank(A) ^ min(m, n) shows that m = n. The system deter¬ 
mining 77^+3 E is therefore a square invertible system. □ 

Lemma 5.7. For all v g and E e i7(curl,7f), 

gradin^l^ v) e curl(iT-l| E) e B^^K). 

Proof. Let e = grad(77^^^^i;) = grad(77Qi;o) where 77 q1;o ^ is as in (5.4). Since 

77q vq is constant along edges of 7f, e must satisfy t • e = 0 along the edges. Moreover, 

(5.18) meK{nlvo) = 0 

due to (5.3b). Hence, integrating by parts on any face / of dK^ 



Summing over all faces / of dK and using (5.18), we conclude that 

Ct • X = 0. 

JdK 

Therefore, to finish proving that e g it only remains to show that (/)n-curie = 0 

for all (/) G P^f^{dK). But this is obvious from the fact that e is a gradient. 

Next, we need to show that a = curl(77^1J.3 E') is in ^^^3(7^). Since it is obvious that 
a G i?p+3(7f), it suffices to prove that 

(5.19) f (/)n-curl(77"!^^3^7;) = 0 

JdK 

for all (f in P^2(^^)* Note that P^2(^^) be orthogonally decomposed into its subspace 
Po{dK) n P^2(^^) its L^(d7f)-orthogonal complement. The latter is a subspace of 
Ppf^{dK) where (5.19) holds (since g Hence it only remains to prove 

that (5.19) holds for 0 in Po{dK) nP^2(^^)- But Stokes theorem shows that (5.19) actually 
holds for all 0 g Po{dK) - cf. (5.10). □ 
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Lemma 5.8. For all v e E e K)^ and a e H^div^ K), 

(5.20a) grad v = ilp+g grad v, 

(5.20b) curl 77™"^ E = curl E, 

(5.20c) div ilpN ^ ~ -E>+2 div a. 

Proof. By Lemma 5.7, 5i = grad(iT®((^^ u) — iTp^ggradu is in Bp™l{K). We will now show 
that 6i satisfies (5.8). Let (j) e Pp{K)^ and consider 

{(p, Si)k = {(p, grad{n^'^^ v - v))k - {(p, n^+sisradv) - gradu ) i ^. 

By Lemma 5.6, 7T^!J.3gradp satisfies (5.2c)-(5.2d), so the last term above vanishes. Inte¬ 
grating the remaining term on the right-hand side by parts, and using (5.2a)-(5.2b), we find 
that 

(5.21) {<p, grad(iTS);.f v - v))k = 0 V</. e Pp{Kf. 

This proves that (0,5 i)r: = 0, i.e., (5.8a) holds. Next, for any ji = F6Pp+i(iL)3, 

we have 


</i,n X (ii> = {F X n,grad(i7®+^ v - v)} - {^i,n x (i7™^^(gradu) - gradu)>. 

The last term vanishes due to (5.2d). Moreover, 

{F X n, grad(7T|^^3^ v - v)) = -(curl F, grad(7T|^^3^ v - v))k = 0 

due to (5.21), so we have proven that (5.8b) holds as well. Hence by Lemma 5.5, 6 i = 0. 
This proves (5.20a). 

To prove (5.20b), we proceed similarly and show that S 2 = curl Ff^^^E — Ffp^^^cmlE is 
zero. By Lemma 5.7, we know that 62 e ^^^3(7^), so if we prove that 

(5.22a) {S 2 ,<P)K = d V</.ePp+l(7^)^ 

(5.22b) (n ■ S 2 , ijfdK = 0 e Pp+ 2 {K). 

then Lemma 5.4 would yield S 2 = 0. To prove (5.22b), 

<n • 52, /^) 3 K = (n ■ (curl(i7“3 E - E) + {I - curlP), 

= <n • curl(iT^"^^ E - E),h) 3 k by (5.7b) 

= <(77^3 E - E) xn, grad^ ijfdK = 0 by (5.2d). 

To prove (5.22a), 

(^2, (P)k = (curl(i7-3 E-E) + {I- il^^g) curl E, (P)k 
= {cml{n;'Pf^E-E),<p)K 

= E-E, curl ^)k + E-E)xn,cP) = 0 


This finishes the proof of (5.22) and hence (5.20b) follows. 


by (5.2e) 
by (5.2c)-(5.2d). 
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Finally, to prove (5.20c), let ^3 = divllp^^a — 77^+2 diver in Pp-^ 2 {K). For any w e 
Pp+ 2 ( 7 F), integrating by parts, 

(^ 3 , w)k = ( div ( i 7^+3 a - u),w)k 

= a -a, grad w)k + (n ■ (ilp^g a - cr),w), 

which vanishes by (5.2e)-(5.2f). Hence ^3 = 0 and (5.20c) is proved. □ 


Proof of Theorem 5.1. The lemmas of this section prove all statements of Theorem 5.1 ex¬ 
cept (5.1b). To prove (5.1b), we use a scaling argument and the commutativity properties 
of Lemma 5.8. Let K denote the unit tetrahedron and let the i7(curl)-Fortin operator on 
iF, defined as above, be denoted by 77^+3 . By the unisolvency result of Lemma 5.6, the 
77(curl, X)-boundedness of the sesquilinear forms and x niTlE) (a con¬ 

sequence of Lemma 2.2), and by finite dimensionality, there is a Co > 0 such that 

for all E E 77(curl,7F). Let Sk : K ^ K he the one-to-one affine map that maps K onto 
a general tetrahedron K. For E : K ^ M^, define ^{E) = {S'^y{E o Sk)^ and HF’Hxcuri ^ 
7^^ IIF’II ^2 + II d'he elementary proofs of the following assertions (i)-(iii) are 

left to the reader. 

(i) E E if and only if $(F) e 

(ii) There are constants Ci,C 2 depending only on the shape regularity of K (but not on 

hx) such that ^ ^ <^ 21 |^III, curl- 

(hi) = 

These three statements imply that 

ClIliT-'-lElllcurl ^ ^ ^ ^06-2||£;|||,„rl- 

While this immediately gives the needed estimate for the 7^-part, namely 
(5-23) II 77^+3 F||^ 2 (^) ^ C||F||^(curi,K), 

we need to improve the estimate on the curl to finish the proof: For this, we use the commu¬ 
tativity property 

24 ) II ^P+3 = II -^p+3 curl £^||l 2 (;^) by Lemma 5.8, 

^ CII curlL;||j^(div,A') by Lemma 5.4. 

The required estimate (5.1b) follows from (5.23) and (5.24). □ 


Before concluding this section, let us illustrate how to use Theorem 5.1 for error analysis 
of DPG methods by an example. 
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Example 5.9 (Primal DPG method for the Dirichlet problem). Consider the broken variational 
problem of Example 3.4: Find g X = H^{Q) x such that (4.1) holds with 

b{{u,an),y) = (grad u, grad y)h + (dn,y}h, Y = 

We want to analyze the DPG method given by (4.3) with 

Xh = {{Wh,n- Th)e X Wh\K e Pp+i{K) and n • Th\eK e tr^i^Rp+i{K) 

for all (tetrahedral mesh elements) X g 
Yh = {vh e H\Qh) : yhW e Pp+3(i^) for all K e Qu]- 

We have already shown in Example 3.4 that the inf-sup condition required for application of 
Theorem 4.1 holds (and Z = {0}). Hence to obtain optimal error estimates from Theorem 4.1, 
it suffices to verify Assumption 3. We claim that Assumption 3 holds with 77 = . Indeed, 

b{ {wh, n-fh),y- n ^+2 y) = (grad Wh, grad(y - y) )h + (n- fh, y - y)h 

= -{Awh, y - y)h + (n-Th- y - PpS y)h = 0 

by applying (5.2a)-(5.2b) element by element. Note that here we have used the fact that 
the discrete spaces have been set so that —/S.Wh\K ^ Pp-i{K) and n • Th — n • gradre/^ is a 
polynomial of degree at most p on each face of dK (i.e., it is in tr^^ i?p+i(7f)), allowing us to 
apply (5.2a)-(5.2b). Applying Theorem 4.1, we recover the error estimates for this method, 
originally proved in [18]. /// 


6. Maxwell equations 

In this section, we combine the various tools developed in the previous sections to analyze 
the DPG method for a model problem in time-harmonic electromagnetic wave propagation. 

6.1. The cavity problem. Consider a cavity 17, an open bounded connected and con¬ 
tractible domain in M^, shielded from its complement by a perfect electric conductor through¬ 
out its boundary dl7. If all time variations are harmonic of frequency cj > 0, then Maxwell 
equations within the cavity reduce to these: 

(6.1a) —lUJiiH -h curlE = 0 in 17, 

(6.1b) —loueE — curl 77 = — J in 17, 

(6.1c) n X E = {) on dl7. 

The functions E^H^J:f2^C^ represent electric field, magnetic field, and imposed current, 
respectively, and i denotes the imaginary unit. For simplicity we assume that the electro¬ 
magnetic properties e and /i are positive and constant on each element of the tetrahedral 
mesh 17/^. The number cj > 0 denotes a fixed wavenumber. In this section we develop and 
analyze a DPG method for (6.1). 

Eliminating 77 from (6.1a) and (6.1b), we obtain the following second order (non-elliptic) 
equation 


(6.2) 


curl ji ^ curl E — u?eE = 
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where / = lujj. The standard variational formulation for this problem is obtained by multi¬ 
plying (6.2) by a test function F g iJ(curl, f2), integrating by parts and using the boundary 
condition (6.1c): Find E g i7(curl, 12) satisfying 

(6.3) cmlE,cm:lF)f 2 -u‘^{eE,F)n = (f,F} 

for any given / g iJ(curl, 12)'. It is well-known [29] that (6.3) has a unique solution for 
every / g iJ(curl, 12)' whenever ou is not in the countably infinite set F of resonances of the 
cavity 12. Throughout this section, we assume uj ^ F. This wellposedness result provides an 
accompanying stability estimate, namely there is a constant C^j > ^ such that 

(6-4) II-£^11 if (curl,r?) ^ C'o; II/II if (curl, 17)' 

for any / G i7(curl, 12)' and E g iJ(curl, 12) satisfying (6.3). Note that the stability constant 
Cuj may blow up as uj approaches a resonance. We continue to use C to denote a generic 
mesh-independent constant, which in this section may depend on CJ,//, and e as well. 

6.2. Primal DPG method for the cavity problem. The primal DPG method for the 
cavity problem is obtained by breaking (6.3). Multiply (6.2) by a (broken) test function 
F G iJ(curl, 12/i) and integrate by parts, element by element, to get 

curl E, curl F)h+(n x ij,~^ curl E,F)h - co‘^{sE,F)h = if,F)h. 

Now set = nx H = x iji~^ curlF^ to be an independent interface unknown which is 

to be found in i7“^/^(div, 12/^). This leads to the variational problem (3.4) with the following 
spaces and forms: 

(6.5a) Xo = i7(curl, 12), Y = i7(curl, 12/^), 

(6.5b) X = il“^/2(div, 512/,), To = i?(curl, 12), 

(6.5c) h^^E.F) = (/i“^ curlF;, curl F)/,-cj^(gF;,F)/,, 

(6.5d) h{H^,F) = iuj{H^,F)h. 

This is the primal DPG formulation for the Maxwell cavity problem. 

The numerical method discretizes the above variational problem using subspaces X/, ^ 
X = Xo X X and defined by 

(6.6a) Xh = {{Eh, nxHh)e i/(curl, f2) x dQh) : n x e Pp+i{Kf, 

and EhW e Pp{Kf for all K e 17/,}, 

(6.6b) Yh = {Eh e i7(curl, Qn) ■ FhW ^ Np+^{K) for all K e Qh}- 

We have the following error bound for the numerical solution in terms of the mesh size 
h = maxxej 7 ^ hx and polynomial degree p ^ 1. 

Corollary 6.1. Suppose {Eh^n x F/,) g X/, is the DPG solution given by (4.3) with forms 
and spaees set by (6.5) and let {E^n x H) e X be the exaet solution of (4.1). Then, there 
exists a C depending only on uj, p, and the shape regularity of the mesh sueh that 

11^ - II if (curl,i?) + 11^ X {H - Hh)\^-i/2(^^x,dQh) 

^ Ch^ (I^liip+i(i7) + I curlF|^p+i(^) -h I curliF|^p+i(^)) . 



26 


C. CARSTENSEN, L. DEMKOWICZ, AND J. GOPALAKRISHNAN 


Proof. To apply Theorem 4.1, we must verify the inf-sup condition (4.5) for the broken form. 
As in the previous examples, as a first step, we verify the inf-sup condition for the unbroken 
form stated in Assumption 1. Given any E e iJ(curl, f7), let /e g iJ(curl, be defined by 
= (/i“^ curl curl — (jj‘^ {eE ^ F) q for all F e i7(curl,f7). Then, (6.2) and (6.4) 

imply 




Cu; sup 

Feif(curl,G) 


Me,f)\ ^ 

II-^11 if (curl,i?) 


i.e.. Assumption 1 holds with cq = C~^. Assumption 2, with c = (jJ~^ is immediately verified 
by (2.11c) and (2.10c) of Theorem 2.3. Hence Theorem 3.1 verifies (4.5) and also shows that 
Z = {0}. The only remaining condition to verify before applying Theorem 4.1 is Assump¬ 
tion 3, which immediately follows by the choice of spaces and Theorem 5.1. 

Applying Theorem 4.1, we find that 


11^ ^^Ilii(curl,i7) + 11^ ^^)II^-V2(div,ai7/,) 


^ C inf 

(G/j, ,71X R}i)eX}i 


11^ ^^Ilif(curl,i7) + 11^ X ^ ^ X ||^-i/2(div,di7/,) * 


Now, H = (icj/i)”^ curl E is an extension to 12 of the exact interface solution nxH. Moreover, 
the interface function n x appearing above can be extended into Xq ^ = {r e i7(curl, 12) : 

Ia: ^ PpEi{F)^}> Since the interface norm is the minimum over all extensions, by standard 
approximation estimates (see e.g., [21, Theorem 8.1]), 


11^ ^^Ilif(curl,i7) 11^ ^ ^^)llif-V2(div,ai7/,) 


^ c 


11^ ^^Illf(curl,i7) + 11^ ^^Ilfr(curl,i7) 


KeQh ■- 


h 


2(51+1) |^|2 


K 


Hn + ^{K) 


+ ++|curlE| 


ijsi+i(K) 


+ 




{K) 


where l/2<5i^p and l/2<52^p+l. Hence the corollary follows. 


□ 


Remark 6.2. Unlike the standard finite element method, for the DPG method, there is no 
need for h to be “sufficiently small” to assert the convergence estimate of Corollary 6.1. This 
property has been called absolute stability by some authors and other methods possessing 
this property are known [23]. 


6.3. Alternative formulations of the same problem. In Example 3.5, we saw that 
a single diffusion-convection-reaction equation admits various different formulations. The 
situation is similar with Maxwell equations. First, let us write (6.1) in operator form using 
an operator A (analogous to the one in (3.10), but now) defined by 

lUJIiH — curlE 
KjOcE -h curl H 


H 


lUja 

— curl 

H 


E 


curl 

luje 

E 
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as A{H^E) = (0, J) for some given J in . However, we will not restrict to right-hand 

sides of this form as we will need to allow the most general data possible in the ensuing 
wellposedness studies. 

We view A as an unbounded closed operator on whose domain is 

dom(A) = G iJ(curl, : n x E = 0 on df?}. 

It is easy to show that its adjoint (in the sense of closed operators) is the closed operator A* 
given by 

—lUjfiH -h cni\E 
—lUjeE — curliJ J ’ 

whose domain is the following subspace of 

dom(A*) = {{H, E) G i/(curl, Qf : E^ = {) on dQ). 

Classical arguments show that both A and A* are injective. To facilitate comparison, we list 
all our formulations at once, including the already studied primal form. 

Strong form: Let x = (iL, E) be a group variable. Set 

Xo = H(cutI, C) X iL(curl, C), T = To = 

Note that Xq = {dom(X),|| • ||if(curi,i 7)}5 dom(X) considered as a subspace of 
iL(curl, (rather than as a subspace of L^(C)^). The Maxwell problem is to find 
X G Xq, given f eYq^ such that Ax = /. This fits into our variational framework (3.2) 
by setting bo to 

bi{x,y) = {Ax,y)Q. 

Primal form for E: This is the same as in (6.5), i.e., with the spaces as set there, 
with b set to b^{H^^E) = E)h and bo set to 

bo{E,E) = cnT\E,cnTlE)a - u‘^{eE,E)Q, 

the electric primal formulation is (3.2) and its broken version is (3.4). 

Primal form for H: Eliminating E from (6.1), we obtain curlcurl iL — = 

curle“^J and a (possibly nonhomogeneous) boundary condition on n x e“^curliL. 
With this in place of (6.3) as the starting point and repeating the derivation that led 
to (6.5), we obtain the following magnetic primal form. Set 

Xo = i/(cnrl,C), To = Xq, 

X = Qh), Y = Hicuil, Hh), 

b^{H,F) = {e-^cm\H,cmlF)h-Lo\yH,F)h, b^{E^,F) = 

With b set to b^ and bo set to 6^, the magnetic primal formulation is (3.2) and its 
broken version is (3.4). 


H 


—lufi curl 

H 


E 


— curl —luoe 

E 
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Ultraweak form: This form is obtained by integrating by parts all equations of the 
strong form. Using group variables x = y = S) and x = set 

(6.7a) Xo = Yq = i7(curl, n) x i7(curl, Q), 

(6.7b) Y = 77(curl, X = dQh) x dQh), 

(6.7c) bo{x,y) = {x,A*y)h, {x,y) = (H-r,n x S)h - (E-r,n x R}h 

and consider formulations (3.2) and (3.4) with bo = b^ and b = W. Note that in 
the definition of 6^, the operator A* is applied element by element, per our tacit 
conventions when using the (•, •)/i-notation. 

Dual Mixed form: Among the two equations in the strong form, if one weakly imposes 
(by integrating by parts) the first equation and strongly imposes the second, then we 
get the following dual mixed form. Set 

Xo = 77(curl, Q) x Yq = Xq, 

X = i7-i/2(curl, Y = H{cm\, ^^h) x 

and consider (3.2) with bo = b^, 

b^{{H,E),{R,S)) = {tojfiH,R)o-{E,cvirlR)h 
+ {lueE + curl H^S)h- 

Its broken version is (3.4) with b set to 6 ^(^t, (i?, S )) = {E^^n x R)h- 
Mixed form: Reversing the roles above and weakly imposing the second equation while 
strongly imposing the first, we get another mixed formulation. Set 

Xo = X H{cutI, 17), Yq = Xq, 

X = i7-i/2(curl, d^h), Y = L^{Qf x 77(curl, 17,,) 

and consider (3.2) with bo = 

b^{{H,E),{R,S)) = {luiiH -curlE,R)h 

+ {zujeE, S)q + {H, curl S)h- 

Its broken version is (3.4) with b set to b^{H-^, iR,S )) = x S}h- 

These form a total of six unbroken and five broken formulations, counting the already dis¬ 
cussed broken and unbroken electric primal formulation. To analyze the remaining formu¬ 
lations, let us begin by verifying Assumption 1 for all the unbroken formulations. To this 
end, label the statement of Assumption 1 with bo set to the above-defined 6 q as “(/)” for 
all I e i7, 5, U, D, M}. Then (analogous to the equivalences in Figure 1 for the elliptic 
example) we now have equivalence of statements (D), (E )^... as proved next. 
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Theorem 6.3. The following implications hold: 


(S) 


(H) 


(E) 


^{D) 


X(^) 


(M) 


Proof. We begin with the most substantial of all the implications, which allows us to go from 
the strongest to the weakest formulation. When there can be no confusion, let us abbreviate 
Cartesian products of L?{Q) as simply L and write W for i7(curl, C) x iJ(curl, C). Clearly, W 
is complete in the i7(curl, C)^-norm. It is easy to see that the graph norms (||x||| + 
and (||x||| + ||74 *x|||)^/^ are both equivalent to the i7(curl, C)^-norm, so W is a Hilbert space 
in any of these norms. These norm equivalences show that the inf-sup condition {S) holds if 
and only if 

(6.8) C'IIxIIl ^ II^^IIl Vx g dom(A). 

{S) (U): The bound (6.8) implied by {S) shows that the range of A is closed. By the 

closed range theorem for closed operators, range of A* is closed. Since A* is also injective, it 
follows that 


(6.9) 


C\\y\\L ^ ||-4*y||L Vy e dom(A*) 


holds with the same constant as in (6.8). This in turn implies that the following inf-sup 
condition holds: 


C'llyllm ^ sup 

xeL 


\{x,A*y)(2 


^yeW. 


Thus, to complete the proof of ([/), it suffices to show that 

( 6 . 10 ) 


i„f sup - inf sup 

xeLy^w fI|l||2/||w xel fIIl||2/||w 


For completeness, we now describe the standard argument that shows that one may reverse 
the order of inf and sup to prove (6.10). Viewing A* : VF ^ L as a bounded linear operator, we 
know that it is a bijection because of (6.8) and (6.9). Hence (A*)“^ : L ^ W is bounded. The 
right-hand side of (6.10) equals its operator norm ||(A*)“^||. The left hand side of (6.10) equals 
the operator norm of the dual of (A*)”^ (considered as the dual operator of a continuous 
linear operator with L as the pivot space identified to be the same as its dual space). The 
norms of a continuous linear operator and its dual are equal, so (6.10) follows. 

(U) (D): Let x = {H, E) and y = (i?, S) be in W = i/(curl, Q) x i?(curl, C). Clearly, 

W is contained in = i7(curl, Q) x . Because of the extra regularity of 5, we may 

integrate by parts the last term in the definition of 6^((i7, E'), (i?, S)) to get that 6^(x, y) = 
h^{x^y). Hence using (?7), 


sup 


K{x,y)\ 

llyllyo 


^ sup 

yEW 


K{x,y)\ 

\\y\\w 


sup 

yEW 


\\y\\w 


^C\\x\\l. 


( 6 . 11 ) 
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Thus, to finish the proof of (£)), we only need to control cmt\H using the last term of . 


IHI, = _ \hE((H,E),(0,S))-(tu^.E,S)\ 

I euii n 11/^ sup ,, ,, sup ,, ,, 

Sel PIIl sel PIIl 

^ sup + \\iujeE\\L. 

yeY,^ 


Using (6.11) to bound the last term, the proof of (D) is finished. 

(D) (H): For any H g i7(curl, f7), set Ih{R) = curl i7, curl i?)^ — R)q. 

We need to prove (i7), which is equivalent to 


(6-12) ||^||if(curl,G) ^ ||^(curl,i7)'- 

Introducing a new variable E = —{i(joe)~^ curliJ, we find that {i(jj)~^iH{R) = —{E^ curli?)^ + 
R)f 2 . Hence 

E), (i?, S)) = {iuj)-^Ih{R) V(i?, S) g i/(curl, Q) x 
Hence (6.12) immediately follows from {D). 

(H) (S): To prove the inf-sup condition (5), it is enough to prove (6.8) for all x = 

(i7, E) G W. Given any E^G e L^(f7)^, the equation Ax = (F, G) is the same as the system 


(6.13a) icjfiH — curlF = F, 

(6.13b) luueE -h curlF = G. 

We multiply (6.13b) by the conjugate of e“^curlF, for some R e F(curl, f7) and integrate 
by parts, while we multiply (6.13a) by —luoR and solely integrate. The result is 

(e“^ curlF, curlF)i 7 + {luo curlF, R)q = (G, e~^ cmi\R)q^ 

-{u‘^IiH,R)q - (zcjcurlF, F)i 7 = -{E,iuR)q. 


Adding the above two equations together, we get the primal form h^{H^R) = £{R) where 
i{R) = (G,g“^ curlF)i 7 — {E^iujR)q. Hence the given inf-sup condition {H) implies 


C'||F||^(curl,i 7 ) 


^ sup 
Reif(curl,i7) 


\b^{H,R)\ 

1111 if (curl, i?) 


Reif(curl,G) ||^||if(curl,G) 


Since \£{R)\ ^ G(||F||i, -h ||G||L)||F||^(curi,i 7)5 fhis provides the required bound for ||F||i,. 
Since curlF is also bounded, equation (6.13b) yields a bound for ||F||i,. Combining these 
bounds, (6.8) follows. 

To conclude the proof of the theorem, we note that the proofs of the implications {U) 

(M), (M) (5), (F) {S) are similar to the proofs of {U) (F), (F) => (F), 

and (F) (5), respectively. □ 


Theorem 6.3 verifies Assumption 1 for all the formulations because we know from (6.4) 
that (F) holds. Assumption 2 can be easily verified for all the broken formulations using 
Theorem 2.3. Assumption 3 can be verified using Theorem 5.1. Hence convergence rate 
estimates like in Corollary 6.1 can be derived for each of the broken formulations. We omit 
the repetitive details. 
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7. Numerical studies 

In this section, we present some numerical studies focusing on the Maxwell example. Nu¬ 
merical results for other examples, including the diffusion-convection-reaction example, can 
be found elsewhere [10, 14]. The numerical studies are not aimed at verifying the already 
proved convergence results, but rather at investigations of the performance of the DPG 
method beyond the limited range of applicability permitted by the theorems. All numerical 
examples presented in this section have been obtained with hp3d^ a 3D finite element code 
supporting anisotropic h and p refinements and solution of multi-physics problems involv¬ 
ing variables discretized compatibly with the i7^(f2)-i7(curl, f2)-i7(div, f2) exact sequence of 
spaces. The code has recently been equipped with a complete family of orientation embedded 
shape functions for elements of many shapes [24]. The remainder of this section is divided 
into results from two numerical examples. 

Example 7.1 (Smooth solution). We numerically solve the time-harmonic Maxwell equations 
setting material data to 


e = /i = 1, cj = 1, 

and Q to the unit cube. To obtain the unit cube was partitioned first into five tetrahedra: 
four similar ones adjacent to the faces of the cube, and a fifth inside of the cube. We have 
used the refinement strategy of [28] to generate a sequence of successive uniform refinements. 
On these meshes, consider the primal DPG method for £', described by (6.5), with data set 
so that the exact solution is the following smooth function. 

El = sinTTTi sin7rT2 sinTrxs, E2 = E^ = {) . 

Instead of the pair of discrete spaces (6.6) that we know is guaranteed to work by our 
theoretical results, we experiment with these discrete spaces: 

(7.1a) Xh = {{Eh,n x Hu) e .ff(curl, Q) x dQu) ■ n x HhU e Np{K), 

and Ehlx ^ Np{K) for all K e 

(7.1b) = {Eh e i7(cnrl, f2^) : e Np^2{K) for all K e 

The observed rates of convergence of the error \\E — Eh\\H{cuT\,f2h) residual norm p 

are shown in Figure 2a. The rates are optimal. This suggests that the results of Corollary 6.1 
may hold with other choices of spaces. Results analogous to those in [3, 8] that allow one 
to reduce the degree of the test space while maintaining optimal convergence rates are 
currently not known for the the Maxwell problem. 

We also present similar results obtained using cubic meshes using i7(curl, f7)-conforming 
Nedelec hexahedron of the first type. Namely, and are set by (7.1) after revising Np{K) 
to Qp-i^p^p{K) X Qp^p_i^p{K) X Qp^p^p_i{K) where Qi^m,n{K) denotes the set of polynomials 
of degree at most and n in the ti,T 2 and directions, respectively. The convergence 
rates reported in Figure 2b are again optimal. 
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Figure 2. Rates for the case of smooth solution and primal formulation. 


Before concluding this example, we also report convergence rates obtained from the ultra- 
weak formulation of (6.7). The discrete spaces are now set by 

(7.2a) Xh = {{E, H, A, A) e x x dQh) x 77“^/^(curl, dQh) ■ 

E\k,H\k e Pp-i{K)^, E^ldK, H^ldK ^ ^pi^) foi" ^ ^ 

(7.2b) Yh = {{F, G) e 77(curl, x Qh) ■ F\k, G\k e Np+2{K) for all K e 
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Figure 3. Rates for the case of smooth solution and ultraweak formulation. 


Recall that the DPG computations require a specification of the y-norm. Using the obser¬ 
vation (made in the proof of Theorem 6.3) that the adjoint graph norm is equivalent to the 
natural norm in i7(curl, 17)^, we set 

||(E,ff)||?.. 2 (ll -®llL2(i7) + Il-f^lli2(fi) + - curl£;||i2(^) + \\iujeE + 

KeQh 

in all computations involving the ultraweak formulation. The results reported in Figure 3 
again show optimal convergence rates. Note that only the errors in the interior variables 
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Figure 4. Construction of Fichera oven. 

E and H (in L^(C)-norni) are reported in the figure. To compute errors in the interface 
variables, we must compute approximations to fractional norms carefully (see [7] for such 
computations in two dimensions). Since the code does not yet have this capability in three 
dimensions, we have not reported the errors in interface variables. /// 

Example 7.2 (Singular solution). To illustrate adaptive possibilities of DPG method and the 
difference between different variational formulations, we now present results from a “Fichera 
oven” problem. We start with with the standard domain with a Fichera corner obtained 
by refining a cube (0,2)^ into eight congruent cubes and removing one of them. We then 
attach an infinite waveguide to the top of the oven and truncate it at a unit distance from 
the Fichera corner, as shown in Figure 4. Setting 

e = /i = 1, cj = 5, 

we drive the problem with the first propagating waveguide mode, 

El = sin7rx2, E2 = E^ = 0 

which is used for non-homogenous electric boundary condition n x E = n x E^ across the 
waveguide section. Analogous to a microwave oven model, we set the homogeneous perfect 
electric boundary condition n x E = 0 everywhere else on the boundary. The above material 
data correspond to about 0.8 wavelengths per unit domain. In all the reported computations, 
we start with a uniform mesh of eight quadratic elements that clearly does not even meet 
the Nyquist criterion. We expect the solution to develop strong singularities at the reentrant 
corner and edges, but we do not know the exact solution. 

First, we report the results from the electric primal formulation, choosing spaces again 
as in (7.1). Figure 5 presents the evolution of the mesh along with the corresponding real 
part of the first component of electric field Ei. Since we do not have the exact solution for 
this problem, we display convergence history using a plot of the evolution of the computed 
residual r] in Figure 6a. (Recall that theoretical guidance on the similarity of behaviors of 
error estimator r] and the error is provided by Theorem 4.1.) Clearly, the figure shows the 
residual is being driven to zero during the adaptive iteration. 
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Figure 5. Every other iterate (iterates 1,3,5,7,9 and 11) from adaptive algo¬ 
rithm applied to solve the Fichera oven problem with the primal formulation. 
Meshes (left) and the corresponding real part of Ei are shown. 
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Degrees of freedom Degrees of freedom 

(a) Adaptivity for primal formulation. (b) Adaptivity for ultraweak formulation. 

Figure 6. Convergence of residual during adaptive iterations for the Fichera oven. 

Next, we solve the same problem using the ultraweak formulation with the spaces set as 
in (7.2) and the y-norm set to the adjoint graph norm as in the previous example. The 
convergence history of the residual norm rj is displayed in Figure 6b. The evolution of the 
mesh along with the real part of Ei is illustrated in Figure 7. 

It is illustrative to visualize the difference between the two different DPG formulations and 
the accompanying convergence in different norms. Figure 8 presents a side-by-side comparison 
of the real part of the electric field component Ei obtained using the primal (left) and 
ultraweak (right) formulations. The same color scale (min = —1, max = 1) is applied to both 
solutions in this figure (whereas the scales of Figures 5 and 7 are not identical). Obviously, the 
meshes are different, but they are of comparable size, so we believe the comparison is fair. The 
primal method, which delivers solution converging in the stronger i7(curl, C)-norm, “grows” 
the unknown solution slower, whereas the ultraweak formulation converging in the weaker 
L^(C)-norm seems to capture the same solution features faster. Both methods ultimately 
approximate the same solution but at different speeds and the ultraweak formulation seems to 
be a winner. Recall that the number of interface unknowns for both formulations is identical, 
but the total number of unknowns for the ultraweak formulation is higher, i.e., the ultraweak 
formulation requires a larger number of local (element-by-element) computations. /// 

8. CONGLUSION 

In addition to presenting the first analysis of DPG methods for Maxwell equations, we have 
presented a technique that considerably simplifies analysis of various DPG methods. The 
idea is to inherit the stability of broken formulations using the known stability of unbroken 
standard formulations and is the content of Theorem 3.1. To obtain discrete stability for the 
Maxwell discretization, a new Fortin operator was constructed in Theorem 5.1. 

We have shown how certain duality identities (proved in Theorem 2.3) can be used to 
verify a critical assumption (Assumption 2) involving an interface inf-sup condition. During 
this process, we have provided a simple technique to prove duality identities (see (2.2c) and 
Lemma 2.2) like 
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Figure 7. Iterates 1,3,5,7, and 9 from adaptive algorithm for Fichera oven 
problem with the ultraweak formulation. Meshes (left) and the corresponding 
real part of Ei are shown. 
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(a) Primal (left) and ultraweak (right) iterate 6. 


(b) Primal (left) and ultraweak (right) iterate 7. 



(c) Primal (left) and ultraweak (right) iterate 8. 



(d) Primal (left) and ultraweak (right) iterate 9. 


Figure 8. Same scale comparison of real part of Ei computed by the primal 
and the ultraweak formulations. A middle slice (through singular vertex) of 
adaptive iterates 6,7,8 and 9 are shown. 


(where the norm on the right hand side is the norm in the space dual to i7“^/^(curl, dK)). 

Finally, the connection between the stability of weak and strong formulations was made 
precise in Theorem 6.3: the wellposedness of one of the displayed six formulations imply the 
wellposedness of all others. Notwithstanding this result, the numerical experiments clearly 
showed the practical differences in convergences among the formulations. 

Before concluding, we mention a few limitations of our analysis and open issues. Conver¬ 
gence results explicit in the polynomial degree p are not obtained by the currently known 
Fortin operators. While construction of local Fortin operators provides one way to prove 
discrete stability, other avenues to reach the same goal (such as the analysis of [27] assuming 
higher regularity, or the analysis of [3] extending the Strang lemma) may prove important. 
Our analysis did not track the dependence on the wavenumber lj. More complex techniques 
are likely to be needed for such parameter tracking [19], including in the unrelated important 
examples of advective singular perturbation problems [10, 12, 14]. 
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